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Abstract 

We show analytically and numerically the existence of double vortex 
solutions in two-Higgs systems. These solutions are generalizations of 
the Nielsen-Olesen vortices and exist for all values of the parameters 
in the Lagrangians considered. We derive analytically the asymptotic 
behavior of the solutions and confirm it numerically by solving the field 
equations. Finally, we show that these solutions can be embedded in 
realistic theories like the two-doublet extension of the standard model. 
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It has recently been shown ||T], |2[ that the Nielsen-Olesen vortex solution P, 
U can be embedded in the standard electroweak model. It was also shown [|] 
that the resulting electroweak vortices are dynamically (but not topologi- 
cally) stable for a finite range of parameters of the electroweak model La- 
grangian. This range of parameters however is not included within the limits 
allowed by current experiments, assuming the simplest form of the Higgs 
sector (one Higgs doublet). 

There are recent theoretical and experimental indications |J that the 
Higgs sector of the standard model may need to be extended to involve at 
least two Higgs multiplets for the electroweak symmetry breaking. This is 
required, for example, in order to achieve consistent unification of the gauge 
couplings on high energy scales 0. The prospect of this possibility raises the 
question of what (if any) soliton-like objects exist in these extended, multiple 
Higgs theories. 

Clearly the Nielsen-Olesen vortex which involves a single Higgs can not 
be directly embedded in multiple Higgs systems unless it is generalized first. 
Recent studies ]7|] have shown the existence of an extra global symmetry in 
realistic two-Higgs doublet models, which breaks during the electroweak sym- 
metry breaking. The breaking of this extra symmetry may lead to topolog- 
ically stable embedded double vortices (the term 'double ' here refers to the 
number of Higgs fields that wind non-trivially and not to the topological 
charge or the flux of the vortices). It is therefore important to study how 
can the Nielsen-Olesen vortex be generalized to multiple Higgs systems, what 
are the properties of the new solutions and how can they be embedded in 
realistic two doublet models. These issues consist the focus of the present 
work. 

Previous studies || attempting to generalize the Nielsen-Olesen solution 
have claimed that such generalized solutions exist only for certain values of 
the parameters of the two-Higgs Lagrangian. These studies, which did not 
make any numerical confirmation of their results, used a very constraining 
ansatz for the asymptotic behavior of their candidate solutions which led to 
incorrect conclusions about the existence of solutions. 

Here, we avoid the use of any ansatz and show that vortex solutions exist 
for any value of parameters in the Lagrangian. We derive the asymptotic 
behavior of these solutions and point out new features that are not present 
in the Nielsen-Olesen vortices. We also obtain the solutions numerically and 
confirm the derived asymptotic behavior. Finally, we consider the realistic 
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two doublet electroweak model and show that the obtained double vortex 
solution can be embedded in this model. 

Consider the Abelian-two-Higgs Lagrangian: 

C = —F^F^ + \\D^,\ 2 + \\D^2\ 2 - V (*i, $2) (1) 

where $1, $ 2 are complex singlets (the Higgs fields), = 8^ — ieA^, F^ v = 
d^A u - d v A v and 

W>i,* 2 ) = x(l^l 2 -^) 2 +T (l$2|2 -" 22)2+ T (l$l|2+l$2|2 -^-" 22)2 (2) 

Consider now the ansatz: 

$1 = ^i/i(r)e ie (3) 
$2 = v 2 f 2 (r)e i6 (4) 

A. = ^ (5) 
er 

Using the ansatz (3)- (5) and a rescaling of the radial coordinate r - 
the field equations of motion obtained from (1) become 

fi(r) + \m - (1 ~ r f r))2 /i(r) - gi(/i(r) 2 - l)A(r) 



/2V) + -£(r) - (1 ^ 2 (r)) / 2 (r) - g 2 (/ 2 (r) 2 - l)/ 2 (r) 



where 



A1 + A3 Ux A2 + A3 V 2 x 

91 = 2 — 92 = ~ 2 — 

e 2 v 2 e 2 v 1 

A3, Wo, A 3/ Vi, 



/viV2e 



-g 3 (/ 2 (r) 2 -l)/ 1 (r) = (6) 

f 2 

-g 3 (/i(r) 2 -l)/ 2 (r) = (7) 

u"(r) - ^V(r) + q 5 (l - ^(r))ii(r) 2 + g 6 (l - v(r))/ 2 (r) 2 = (8) 



Single valuedness for the fields $1, $2 imply the following boundary condi- 
tions for the system (6)- (8) 

fx -> f 2 -> u -> for r -> (9) 

and 

/l -»• 1 h -»• 1 u -»• 1 for r -> 00 (10) 

It is easy to see that for r 1, /i(r) ~ h{r) ~ r and v(r) ~ r 2 . 

The asymptotic behavior of the fields may also be obtained in a straight- 
forward way. Define Sfi, 5f 2 and Sv by 

fi^l + 5fi v^l + Sv for r>l (11) 

where i = 1,2. By keeping only lowest order terms in Sf\, Sf 2 , the system 
(6)- (8) becomes: 

5tf(r) + hf{(r)- 5 -^-2q 1 8f 1 (r)-2q 3 8f 2 (r) = (12) 
SfS (r) + -8m- 6 -^-2q 2 5f 2 (r)-2q 4 5f 1 (r) = (13) 

S v »( r )-~5 v >(r) + (q 5 + q 6 )5v(r) = (14) 
r 

Equation (16) shows that Sv decays exponentially. This implies that the sys- 
tem of (12) and (13) becomes mathematically similar to a system of coupled 
harmonic oscillators driven by the same 'force'. This is a simple problem but 
we will go through it in some detail in order to show that there is always a 
solution and thus resolve the controversy with Ref. ||. Notice that the term 
&v( ^2 has been kept in (12) and (13). Neglecting this term, as was done in 
the original paper of Nielsen-Olesen , leads to an asymptotic behavior for the 
fields, which is incorrect for large values of the self coupling when this term 
dominates. For a proper treatment of this term in the Nielsen-Olesen case 
see Ref. M. 

Clearly the solution to (14) may be approximated for r ^> 1 by 

5v = r T e- ar {c v 1 +c v 2 r- 1 ) (15) 



Using this ansatz in (14) leads to 



This gives the behavior of the 'driving force' in the system (12), (13). The 
homogeneous solution of this system may now be obtained by using an ansatz 
of the form 

Sfi = c{r a e-^ (17) 

which leads to 

c{((3 2 - 2 Ql ) - c f 2 2q 3 = (18) 

-c{2q A + 4((3 2 - 2q 2 ) = (19) 

Demanding that this system has a solution implies that the secular determi- 
nant vanishes which in turn leads to 



The corresponding ratio -j is 



P± = [(<Zi + fc) ± yj(qi ~ g 2 ) 2 + Msqi (20) 

J 



,cL (?2 -qi) ± (q 2 ~ qi) 2 + 4:q 3 q 4 

(-)± = 7T = 7± (21) 

c{ ^3 

By definition q 3 q 4 < qiq 2 (for Ai,A 2 ,A 3 > 0) which implies j3± G 3? and 
7_ < < 7+. Therefore, the homogeneous solution to the system (12), (13) 
is 

5f* = (cf + e _/?+r + c{_e-P- r )r a (22) 
5tf = (c^+e"^ + 4-l-e~ fi - r )r a (23) 

where a is a constant that can be determined numerically. The particular 
solution of the system is of the form 

Sff = d{ e — (24) 
where i — 1,2. The ansatz (24) leads for r 3> 1 to the inhomogeneous system 

d{{Aa 2 - 2qy) - d f 2 2q 3 = {c[) 2 (25) 
-d{2q 4 + d f 2 (4a 2 -2q 2 ) = (c^) 2 (26) 
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where c\ is denned in (15). The constants d{,d 2 can be obtained by solving 
the system (25), (26) which always has a solution. Therefore, the form of 
Sfi, Sv is: 

„-2<rr 

Sf, = (c{ + e- p+r + c{_e" /3 - r )r a + d{ (27) 



r 

Sf 2 = {c{ +1+ e~?+ r + c{_^e^- r )r a + d[— — (28) 

Clearly, the term with the lowest value of the exponent will dominate. We 
may therefore distinguish four cases 

1. Ai 7^ A 2 (or V\ 7^ v 2 ) and /?_ < 2a. In this case the term proportional to 

e -/3_r Joniinates and we have 

5^ = c{_e~ p - r r a (29) 
Sf 2 = c{_7_ e -' 3 - r r a (30) 
§ v = c \e- ar r^ (31) 

Notice that since 7_ < 0, Sf 2 is negative in this case and therefore 
/ 2 (r) is approaching its asymptotic value (=1) from above. This is a 
new feature which does not appear in the Nielsen-Olesen solution for 
any value of the parameters. We have checked the validity of this point 
numerically (see below). 

2. Ai 7^ A 2 (or v\ ^ v 2 ) and (3- > 2a. In this case the term proportional to 

e -2ar dominies and 

f e- 2ar 

Sfi = d (32) 

r 

Sv = c?e _< M (33) 

where i — 1, 2. 

3. Ai = \ 2 , Vi = v 2 and f3 + < 2a. This implies by symmetry that Sfi = Sf 2 

and therefore, since 7_ < 0, we must have c{_ = 0. Therefore 

Sfi = c{ + e-?+ r (34) 
Sv = cle'^A (35) 

since j + — 1, a — \[2 for Ai = A 2 and V\ = v 2 . 
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4. Ai = A2, V\ — V2 and (3 + > 2a. The dominant term now is proportional 
to e~ 2crr and 5f\ = 6/2 by symmetry. Thus 

e -2V2r 

Sfi = d{ (36) 

r 

5v = c\e-^ r r^ (37) 

We have verified the above cases by plotting ln ^ s J 1 ^ and V s r for 

several different values of parameters. In Figures 1 and 2 we show those 
plots for parameter values Ai = A 2 = A 3 = 1 (case 3) and Ai = A 3 = 1, 
A 2 = 2 (case 1) respectively. The asymptotic values for all fields are consistent 
with the corresponding predictions of exponents. Indeed, for Figure 1 the 
relevant exponent for 8f\, 6/2 is (3 + = 2.45 (case 3) while for Figure 2 the 
corresponding exponent is /?_ = 1.7 (case 1). In both cases we had v\ = v% 
which implies that the relevant exponent for 5v is a = \pl. Notice also the 
singularity of ln ^ 5 J 2 ^ in Figure 2, shown as a sharp minimum due to the 
discrete nature of the numerically constructed plot. This is due to the fact 
that / 2 (r) crosses the / 2 (r) = 1 line (leading to 5f 2 = 0) in order to approach 
its asymptotic value / 2 = 1 from above as predicted by (30) (since 7_ < 0). 

It is instructive to study how does the above analysis get modified in 
the case of more general two-Higgs potentials appearing in realistic cases. 
Consider for example the potential: 

F($i, $ 2 ) = V ($i,$ 2 ) + A 4 (|$i| 2 |$ 2 | 2 - |$;$ 2 |) + A 5 |$;$2 - v lV2 \ 2 (38) 

where Vo($i,$ 2 ) is defined in (2). Using the same ansatz as in (3)-(5) we 
obtain the field equations for /i(r), / 2 (r) and v(r) as: 

f"(r) + l -f[(r) - (1 ~ r f r))2 /i(r) - ?i(A(r) 2 - l)A(r) - 

? 3 (/ 2 (r) 2 - l)/x(r) -p(/i(r)/ 2 (r) - l)/ 2 (r) = (39) 

fi(r) + -f a (0 - (1 ~f r))2 / 2 (r) - g 2 (/ 2 (r) 2 - l)/ 2 (r) - 

g 3 (/i(r) 2 - l)/ 2 (r) -p'(/i(r)/ 2 (r) - l)A(r) = (40) 

v "( r ) _ I v '( r ) + q5 (l _ ^(r))/^) 2 + g 6 (l - ^(r))/ 2 (r) 2 = (41) 
r 

where 

P=— — an d P = — — 42 

t»i e 2 v 2 
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and qi as denned previously. It is easy to show that the equations for Sfi, 5v 
may be obtained from (12)-(14) by substituting 



2 Ql -> 2 qi +p 2q 2 ^2q 2 +p' (43) 
2q 3 -> 2q 3 +p 2q 4 -> 2q 4 + p' (44) 

and the analysis proceeds in exactly the same way but with different param- 
eters. 

We now show that the obtained solutions can be embedded in realistic 
theories. The main steps are similar to those in Ref. JIJ which showed that 
the Nielsen-Olesen vortex is a solution of the equations of motion obtained 
from the standard electroweak Lagrangian involving a single Higgs doublet. 
We therefore follow the notation of Ref. |TJ with minor modifications. 

Consider the bosonic sector of the two-Higgs doublet electroweak La- 
grangian: 

L = L W + L B + L* u * a -V($i,$ 2 ) (45) 

with 

L w = ~ l -G^ a G^ a (46) 
L B = -\f b , u F b ^ (47) 

L* li#2 = El^ii 2 ^j^\{d x -\igr a W^-\ig'B x )^\ 2 (48) 
i=i i=i A A 

where G^ va , F BfMU are the usual SU{2) and U(l)y gauge tensors and W x , B x 
are the corresponding gauge fields. Also y($i,$ 2 ) is given by (38) with $i, 
$ 2 being complex doublets. 

The field equations obtained from (45) are 

2 j 

d v G^ a - gs abc G» vb WZ = Yl 7;9{^W a D^i - (D^V^) (49) 

i=i ^ 

d v F^ = Y'-g'^lD^-iD^^) (50) 

w = j_m^A (51) 



7 



which are almost identical to the ones shown in Ref. |T| apart from the 
two-Higgs potential term. Define now e, 9 W , Z^ and A 11 as usual i.e. 



g = e cos 9 W g' = e sin 9 W 
Z» = cos9 w W^ -sm9 w B^ 
A" = sin 9 W W^ + cos 9 W B" 



(52) 
(53) 
(54) 



and consider the ansatz 



Z^ 



yyixl _ iy^2 = q 



(55) 
(56) 



er 




(57) 



(58) 



It is straightforward to show that the ansatz (55)- (58) substituted in (49)- 
(51) leads to the equations (39)-(41) for fi(r), /2O") and v(r). Therefore, 
the double vortex solution is also a solution of the field equations in the two 
Higgs doublet electroweak model. 

It has been argued that a certain type of the embedded double vor- 
tex is topologically stable due to an extra global symmetry that amounts 
to the freedom of independent U(l) phase transformations of the two dou- 
blets. The breaking of this symmetry during the electroweak transition leads 
to stability of the embedded double vortices. In the Lagrangians we have 
considered however, the presence of the A5 term does not allow such global 
transformations and it explicitly breaks the extra global symmetry. Even in 
the absence of this term the investigation of the stability of the embedded 
vortices needs special attention. In fact, instabilities may be present even in 
the topological case towards repulsion of the centers of the two vortices. This 
is similar to the Nielsen-Olesen n-vortex (n being the winding number) which 
even though topological is unstable towards decay to n 1-vortices for large 
self coupling. In the case of double vortices it would cost logarithmically 
infinite energy to place the two centers infinite distance apart. This is due to 
the fact that the gauge field can not effectively screen the angular gradient 
energy of the two centers when they are separated. This however leaves open 
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the possibility of repulsion at small distances which would be favoured by 
the two-Higgs potential term. These issues of stability are currently under 
investigation. 
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Figure Captions 

Figure 1: The dependence of ln ^ x ^ on r for x = 8 fx, x = 5f 2 and x = 5v 
with Ai = A 2 = A 3 = 1, V\ = v 2 (case 3). The predicted asymptotic values 
are f3 + = 2.45 for 8f%, Sf 2 and a = v2 for 5v. Clearly, the numerically 
obtained asymptotic values are in good agreement with the predicted ones. 

Figure 2: Same as Fig. 1 with Ai = A 3 = 1, X 2 = 2, vi = v 2 (case 1). 
The predicted asymptotic values are /3_ = 1.7 for 8fi, 5f 2 and a = \[2 for 
5v. 
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